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^ ■ Abstract 



The previously developed renormalizable perturbative 1 /N-expansion in higher 
^ ' dimensional scalar field theories is extended to gauge theories with fermions. It is 

based on the 1/A^j-expansion and results in a logarithmically divergent perturba- 
tion theory in arbitrary high odd space-time dimension. Due to the self-interaction 
of non-Abelian fields the proposed recipe requires some modification which, how- 
ever, does not change the main results. The new effective coupling is dimensionless 
and is running in accordance with the usual RG equations. The corresponding 
beta function is calculated in the leading order and is nonpolynomial in effective 
coupling. The original dimensionful gauge coupling plays a role of mass and is also 
logarithmically renormalized. Comments on the unitarity of the resulting theory 
are given. 



1 Introduction 

In our previous paper [1] (hereafter [I]) we constructed renormalizable perturbative 1/N 
expansion in scalar field theories in higher dimensions. It was demonstrated that the 
1/N expansioiil] contrary to the usual perturbation theory is renormalizable [HI [6], leads 
to a finite number of logarithmically divergent terms and can be made finite by renor- 
malization of the fields and a new dimensionless coupling. The original dimensionful 

^for review of the 1/N expansion see [21 [31 H] 
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coupling does not serve as an expansion parameter anymore and plays the role of mass 
which is also logarithmically divergent and multiplicatively renormalized. Within the 
dimensional regular izat ion technique [7] we performed the renormalization procedure 
in scalar theories in arbitrary odd space-time dimension and calculated a few terms of 
the 1/N expansion. Even dimensions, in principle, can also be treated by this method; 
however, they lead to some complications due to the appearance of log terms. We did 
not consider them in [I] and do not do it here. 

The main recipe of the 1/N expansion is the following (see for example [8]): one 
considers a theory with an N-component field interacting with some real or auxiliary 
singlet field and calculates the propagator of this singlet field. If one divides the coupling 
of their interaction by V^, then any simple virtual loop in this propagator is suppressed 
by a factor of 1/N which is cancelled by an opposite factor coming from the closed loop 
of the N-component field. Thus, each of this bubbles has no 1/N suppression and in 
the leading order one has to sum them up. The resulting "zeroth" order propagator 
contains the polarization operator, which is an increasing function of momenta, in the 
denominator. This new propagator has to be used now in all Feynman diagrams order 
by order in 1/N. Since it decreases faster than the usual one, it improves the convergence 
of the integrals and leads to logarithmically divergent diagrams even in arbitrary high 
extra dimensions. 

The aim of this paper is to demonstrate how the advocated procedure works in gauge 
theories which contain some peculiarities due to the gauge invariance. We show that 
indeed in the case of an Abelian theory the procedure is straightforward while in the 
non-Abelian case it needs some modification due to the presence of the triple and quartic 
interactions of gauge fields. 

When summing up the vacuum polarization diagrams to the denominator, one faces 
the problem of unitarity due to the presence of the imaginary part in the polarization 
operator and possible poles in the complex plane. This problem is common to any 
realization of the 1/N expansion and deserves a special treatment [6], [9]. We briefly 
comment on this problem below. 

2 QED 

Let us start with the usual QED with Nj fermion fields in D dimensions, where D takes 
an arbitrary odd value. The Lagrangian looks like 

1 1 ---e-- 

C = --{d^A^- d^Af,f - —{d^A^f + ii)idi)i-mi)ii)i + -j=i^^^ (1) 

According to the general strategy, we now have to consider the photon propagator. 
Since due to the gauge invariance the polarization operator is transverse, it is useful to 
consider a transverse (Landau) gauge. This is not necessary but simplifies the calcula- 
tions. Then in the leading order of the 1/N expansion one has the following sequence 
of bubbles (see Fig.l) 
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Figure 1: The chain of diagrams giving a contribution to the A field propagator in 
the zeroth order of the 1/iV/ expansion 



summed up into a geometrical progression. The resulting photon propagator takes the 
form similar to that for an auxiliary a field in [I] 



D,u{p) ^,1^9^" p2 J i + e2jp)(_p2)D/2-2' (2) 
where 

..^._ r\D/2n2-D/2) 

jy^) 2[^/2i-ir(D)7r^/2 

and we put m = for simplicity. 

This practically coincides with the expression obtained in scalar theory and all the 
following steps just repeat those in the latter. As was explained in [I], we change the 
normalization of the gauge field A^ ^^/^ ^^"-^ introduce the dimensionless coupling 
h associated with the triple vertex, so the effective Lagrangian takes the form 
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+iiljidijji - mi/jii/ji + ^"(^ ^IJiA^lJi. (3) 



This new dimensionless coupling h enters into the gauge transformation and plays the 
role of a gauge charge. The old coupling e, on the contrary, is dimensionful and acts as 
a mass parameter in a gauge propagator. Since the coupling constant h is dimensionless 
the effective Lagrangian ([3]) when omitting the first term is conformal as considered 
in [To] where the theory was taken in D = 3. 

Again, one has the modified Feynman rules with the photon propagator that de- 
creases in the Euclidean region like l/{p^)^^'^~^ , thus improving the UV behaviour in a 
theory. The only divergent graphs are those of the fermion propagator and the triple 
vertex. They are both logarithmically divergent for any odd D. The photon propagator 
is genuinely finite and may contain divergencies only in subgraphs. One basically has 
the same graphs as in [I] for a scalar theory but with solid lines being the fermion ones 
and the dashed lines being the photon one. 

The only difference (or simplification) comes from the Furry theorem and the gauge 
invariance. Namely, all triangles with three photon external lines vanish due to the 
Furry theorem and the gauge invariance which connects the fermion propagator with 
the triple vertex implies that Zi = Z2. This relation holds in the expansion like 

in the usual PT. Thus, using the notation of [I], in the leading order one has 

r(D)(D-i)(2-D/2) 
2^/2r(2-D/2)r(Z}/2 + i)r2(z}/2)' , ^ ^- yj 



The same results were obtained in [TT] where the author calculated the anomalous 
dimensions at the d-dimensional critical point where the fields obey asymptotic scaling 
and are conformal. This leads to the following renormalization constants in the leading 
order in l/Nj: 

1 Ah 

and, consequently, hs = h. Hence, in odd-dimensional QED in the leading order of the 
expansion one does not need the coupling constant renormalization; only the wave 
function renormalization remains. This means that the coupling is not running. 

In the second order one again has the same diagrams as in a scalar theory but with 
vanishing triangles. The renormalization constant in the second order is also essentially 
simplified compared to the scalar case and looks like 

I Ah 11 A^h'^ 1 ^ 

As explained in [I] the leading pole term ~ 1/ s^N'j in ([6]) can be evaluated alternatively 
via the pole equations [12] from the 1/eNf one with the same result. 

Like in the scalar case [I] the original dimensionful coupling e is not an expansion 
parameter anymore, but plays a role of a mass and is multiplicatively logarithmically 
renormalized. The leading order diagrams are shown in Fig.2. 

"^^'(^y^ ™<0>™ 
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Figure 2: The first order diagrams giving a contribution to the 1/e^ renormalization 
in the 1/Nf expansion 

They give the following contribution: 

Diag.a =^ rr^-^; Diag.b =^ rr^-^i Diag.c =^ 0, (7) 



eNf{l + hy ' ^' eNf{l + hy' 
r{D + l){D/2 - 1){D - l)\2 - D/2) 0^^0/2-9 

2^/2+ir(2-D/2)r(D/2 + 2)r2(z}/2)' d/2{d/2-i){d -i) ■ 

So one has 



The solution of the RG equation for 1/e^ with fixed h gives the momentum dependence 
of 



e2 el \pIJ 
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with the anomalous dimension 

_ T{D){D - l){D/2 - 2){D - 3){D + 2){D - 6) 
^ ~ 2^/2+ir(D/2 + 2)r2(D/2)r(2 - D/2)Nf (l + /i)3' 

The sign of 7 depends on D. For D = 5, 7 7 > 0, for D = 9 7 < and then alternates 
with every odd D. Notice that eq.Q reminds the power law behaviour of the initial 
coupling in extra dimensions within the Kaluza-Klein approach [13j . 



3 QCD 



3.1 The effective Lagrangian 

Consider now a non-Abelian theory. Here one has some novel features due to the 
presence of the triple and quartic gauge vertices and the ghost fields. Similar to ([1]) we 
write down the Lagrangian for the gauge fields and Nf fermions as 



C 

where 
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2a 



duAl + 
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A" 



Like in QED we choose the Landau gauge and sum up the fermion bubble diagrams into 
the denominator of the gauge field propagator 



G 



ah 



{9 



HP 



p2 l+^2Jp)(_p2)D/2-2^ 



(10) 



where the coefficient f{D) differs from the Abelian case only by the color factor T{R) 

V^{D/2)V{2- D/2) 



2l^/2]-ir(D)7r^/2 



-TiR) 



and again we put m = for simplicity. 

In the non-Abelian case, contrary to the Abelian one one has the triple and quartic 
self- interaction of the gauge fields. These vertices, which are suppressed by 1/ ^/Nj and 
l/Nf, respectively, obtain loop corrections of the same order in 1/Nf. The effective 
vertices in the leading order are given by the diagrams shown in Fig. 3 and 4. 





Figure 3: The diagrams giving a contribution to the A^ term in the zeroth order 
of the 1 / Nf expansion 
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Figure 4: The diagrams giving a contribution to the term in the zeroth order 
of the 1/Nf expansion 

Thus, besides the modification of the gauge propagator one has the modified vertices. 
The effective Lagrangian in the case of vertices is not given by a simple local expression 
due to complexity of the loop diagrams. So we keep it in the form of the diagrams 
which have to be evaluated in integer dimension. Due to the rules of dimensional 
regularization they are finite for any odd D, otherwise one has to redefine them. What 
is crucial, however, is that there are only three diagrams which have to be redefined. 
Hence, after rescaling the gauge field A^/g one obtains the following effective 

Lagrangian: 



C 



eff 



1 

2ag^ 



(11) 
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Notice that dimensionful coupling g drops from all terms except for the first one and is 
not an expansion parameter anymore. 

Calculating the degree of divergence after summing up the diagrams of the zeroth 
order, similar to the scalar case and QED, one has only four types of logarithmically 
divergent diagrams: the fermion and the ghost propagators, the fermion-gauge-vertex 
and ghost-gauge-ghost vertex. The gauge propagator as well as pure gauge vertices are 
finite and may contain only divergent subgraphs. 

The next step is the introduction of a dimensionless coupling h. Here one should be 
accurate since this coupling enters not only into the triple gauge-fermion vertex, but due 
to the gauge invariance should be present in gauge and gauge-ghost vertices. It should 
be the same in all three of them. In the case of a gauge theory, the coupling h enters 
the gauge transformation and acts as a gauge charge of the fermion and gauge fields. 

When constructing the Feynman diagrams, one reproduces the one-loop cycles that 
are already present in the effective Lagrangian f|T2|) but with additional factors h. In 
the scalar or QED case, this happened only for the propagator, but here it is also true 
for the vertices. As a result, the final expression for the effective Lagrangian takes the 
form 



C 



1 




e// = - - (O +X\ + XX ) (1+^) (12) 



h 



2ag^ 



where 



cabc Ab jz 
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3.2 Properties of the l/Nf expansion 

Consider now the leading order calculations. We start with the terms for the 

fermion and the triple fermion-gauge-fermion vertex. The diagrams are shown in Fig. 5. 
The first two are the same as in QED. The third diagram contains new effective vertex 
which includes the usual triple vertex and the fermion triangle. The usual vertex does 
not give a contribution since it is finite by a simple power counting. At the same time, the 
fermion triangle is momentum dependent and the resulting diagram is logarithmically 
divergent. 




a b c 

Figure 5: The leading order diagrams giving a contribution to the if) field propagator 
and the triple vertex in 1 / Nf expansion 



Calculating the singular parts of the diagrams of Fig. 5 in dimensional regularization 
with D' = D — 2e one finds 

Diag.a I ^ A, Diag.b -^FTVT^' Diag.c -^FTVT^' f^^^) 

eNf 1 + h eNf 1 + h eNf 1 + h 

V{D){D-1){2-D/2)Cf Cf-Ca/2 (1 - D/2)Ca 

2D/^V{2- D/2)V{D/2 + l)V\D/2)T' Cp 2{2-D/2)Cf ' 



which is again in agreement with [H]. Notice that the third diagram is proportional to 
h/{l + h) instead of h'^/{l + as in the scalar case. The reason is that now we have 
an effective triple gauge vertex proportional to \/h{l + h) instead of y/h that cancels 
one factor oi h/(l + h). 

Therefore, in the leading order in the expansion the renormalization constants 
take the form 

= 1-^^, (14) 
1 iB + C)h 

2' = i-^TTTftr- (15) 

= zfZ- = 1 _ ^M±£±^. (16) 

To check the gauge invariance, we calculated the renormalization of the coupling 
through the gauge-ghost interaction. The leading diagrams are shown in Fig.6. 
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a 



b 



c 



Figure 6: The leading order diagrams giving a contribution to the ghost field propagator 
and the triple vertex in expansion 



Calculating the singular parts of the diagrams in dimensional regularization one finds 

l/i , 1 h , 1 h , 

D^ag.a ^JTI^ , D^ag.b ^ ^^^5 , D^ag.c ^^^^ , 

..^ rmD-i)CA 

2D/2+ir(2-D/2)r(D/2 + i)r2(D/2)r' u, ^ u, ^ ,) 

which gives the following renormalization constants in the ghost sector 

Zi = 1, (18) 
~ 1 1 A'h 

= Z?Z-=1-— — . (20) 

One can see that the following relation holds: 

A + B + C = A' + B' + C, (21) 

which follows from the gauge invariance. 

We look now at the next-to-leading order to compare it with the scalar case. The 
corresponding diagrams for the fermion propagator are shown in Fig. 7. They require 
some explanation. The first line of diagrams in Fig. 7 is obtained from the one- loop 
diagrams of Fig. 5 by inserting into the vertex or the fermion line of the one-loop divergent 
subgraphs from Fig.5. For example, the diagram d in Fig. 7 is the diagram a from Fig.5 
with divergent one-loop subgraph c from Fig.5 substituted instead of the initial vertex. 
The second line of the diagrams in Fig. 7 is obtained from the "forbidden" diagram of 
Fig. 8 by inserting the same one-loop divergent subgraphs from Fig.5. The diagram e is 
the diagram of Fig. 8 with insertion of the subgraph a from Fig.5 into the fermion line 
(see Fig. 9) and the diagram g comes from the insertion of the subgraph c from Fig.5 
instead of one of the vertices in the fermion loop (see Fig. 10). 
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A 





f 



Figure 7: The second order diagrams giving a contribution to the fermion propagator 
in the expansion 




Figure 8: The "forbidden" loop diagram 




Figure 9: The diagram e from Fig.7 as a result of insertion of the diagram a from Fig. 5 
into the fermion line. 




Figure 10: The diagram g from Fig.7 as a result of insertion of the diagram c from Fig.5 
instead of one of the vertices in the diagram from Fig. 8. 

All these diagrams are double logarithmically divergent, i.e., contain both single and 
double poles in dimensional regularization. We calculated the leading double poles after 
subtraction of the divergent subgraphs, i.e., performed the i?'-operation. The answer is: 



Diag.a 



1 



2u2 



Diag.b 



1 ABh^ 



2u2 



Diag.c 



Diag.d 
Diag.f 



1 ACh^ 
1 2ABh^ 



Diag.e 



Diag.g 



1 2A^h^ 
1 2ACh^ 



(22) 



We performed also the calculation for the fermion-gauge-fermion vertex but do not 
present the diagram-by-diagram result because of the lack of space and give only the 
final answer. 

Having all this in mind we come to the final expressions for the Z factors in the 
second order of the 1/iV/ expansion in the fermion sector: 



1 {B + C)h^ 1 



eNf l + h 



e^Nj \2 



3{B + Cyh^ A{B + C)h^ 



2 {B + Cfh^ 2 A{B + C)h^ \ 1 



(23) 



1 - 



Ah 1 

+ 



7 

3 A^h"^ ^ A{B + C)h^ 



eNf l + h e'^N] \2{l + hy (l + h) 



2 A{A + B + C)h' \ .^.J_ 
3 + J^^eN^ 



The same calculation in the ghost sector gives 
Zi = 1, 



3 A'^h^ 



eNf l + h e^Nf \2 (1 + h) 



A'h 1 

+ 



/ 



2A'{A + B + C)h^ 
3 {1 + hf 



(24) 



(25) 



+ 0(^)(26) 



Notice the absence of the ghost-gauge-ghost vertex renormalization. 

The final second order expression for the coupling renormalization calculated in both 
ways having in mind relation (!?!]) is 



Zh 



^ 1 2{A + B + C)h ^ 1 



eNf 



l + h 



e-N] 



{A + B + CfK' A{A + B + Cfh'^ 



{l + hf 



1 



eN]' 
(27) 

Like in the scalar and QED case, one can also calculate the renormalization of the 
original coupling g"^. The leading order diagrams are shown in Fig. 11 which give the 
following singular parts like in 



^ -<D- -TV 




a 



Figure 11: The first order diagrams giving a contribution to the l/g"^ renormalization 
in 1/A^/ expansion 
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r(Z} + l)(Z}/2-l)(D- 1)^(2 -D/2)C^ D^ + D/2-9 Cf-CaI2 

2D/2+iY{2- D/2)T{D/2 + 2)T^{D/2)T{R)' D/2{D/2 - 1){D - 1) Cp 

A{D/2f - Q{D/2f + 18{D/2y - Q7{D/2f + 85(£'/2)^ -19D + Q Ca 

2{D - 1)2(1 - D/2)2(2 - D/2)D CV ' 

D^-D^/2-2D + l Ca 
D{1-D/2){2-D/2){D-1)^Cf ' 

The corresponding renormalization constant looks like 

7 1 ^ {F + E + G + H)h' 

3.3 Renormalization group in l/Nf expansion 

Having these expressions for the Z factors one can construct the coupling constant 
renormalization and the corresponding RG functions. One has as usual in dimensional 
regularization 

hs = {^.rhzlz,^ = {^^r(h + f:^^^^^^ (30) 



n=l 



z. - 1 + (31) 

n=l 

where the first coefficients a„ and can be deduced from eqs. ( I23|2^25II26I) . 

This allows one to get the anomalous dimensions and the beta function defined as 

7(/.iV/) = V;^logZ = /4c„ (32) 
l3(h,Nf) = 2h('a+l2) = (h^-l)ai. (33) 
With the help of eqs. fl23ll24p one gets in the leading order of the 1/A^/ expansion 

l2{h,Nf) = ^,(h,Nf) = -—^4^^, (34) 

' fJ Nf{l + h)^' ' Nf (l + /i)2 ' ^ ^ 

l2ih,Nf) = -^(YTIF' 7i(/^,iV,) = 0(i^), (35) 
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It is instructive to check the pole equations [12] that express the coefficients of the 
higher order poles in e of the Z factors via the coefficients of a simple pole. For one 
has, according to (!24l) . 



1 Ah 

cM ^ (37) 
1 /3 A^h^ , A{B + CW 2 A{A + B + C)h' \ 



At the same time, the coefficient C2 can be expressed through ci via the pole equations 

as 

Integrating this equation one gets for C2 the expression coinciding with (l38l) which was 
obtained by direct diagram evaluation. 

We have also checked the pole equations for the renormalized coupling. From eq. ( l27|l 
one gets the coefficients of the coupling constant renormalization factor 

"'f"'^^) = -W, i + A • <«> 

1 A(-4 + S + C)'^ft' 4(/l + B + C)';i*\ 
«2(ft,A'/) = ]vj(^ (l + ft)2 3 + j- 

At the same time, from the pole equations one has 

The coefficient 02 evaluated in this way coincides with fHTj) . 

Equation fl27|) gives us the sign of the beta function. In the leading order one has 



^ = ^ih) = nD){D-l)CA ^ 

dt ' 2^/^r{2- D/2)r{D/2 + l)r\D/2)NfT{l + hy' ^ ' 

which means that (3{h) < for D = 5, j3{h) > for D = 7 and then alternates with D. 

Notice that in QCD, contrary to QED, all Feynman diagrams contain group factors 
so that the actual expansion parameter becomes N^/Nf, thus requiring that this ratio 
is small. Of course, in non-Abelian theories the 1/Nc expansion would be preferable, 
since it accumulates the interactions of the gauge fields, however, in this case already 
the lowest approximation consists of all planar diagrams and is not known [15] . 



4 Conclusion 

We conclude that in higher dimensional gauge theories like in the scalar case despite 
formal non-renormalizability it is possible to construct renormalizable 1/A^/ expansion 
which obeys all the rules of a usual perturbation theory. The expansion parameter is 
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dimensionless, the coupling is running logarithmically, all divergencies are absorbed into 
the renormalization of the wave function and the coupling. The original dimensionful 
coupling plays a role of mass and is renormalized multiplicatively. Expansion over this 
coupling is singular and creates the usual nonrenormalizable terms. 

Properties of the expansion do not depend on the space-time dimension if it is 
odd. In even dimension our formulas after subtraction contain a logarithm which creates 
some technical problems in calculations. We plan to consider these theories later. 

There is one essential point that we omitted in our discussion, namely, the analytical 
properties of the gauge field propagator and the unitarity of a resulting theory. As one 
can see from eqs.(l2]) and ( |T0|) . the propagator of the gauge field contains a cut starting 
from = in the massless case or from 4m^ in the massive one. It may also contain 
poles in complex momentum plane. This imposes the question about unitarity of a 
resulting theory. 

The unitarity of such an expansion was shown in [8] at the tree level. In the case 
of loops, one can show that all the cuts imposed on diagrams when applying Cutkosky 
rules [16] in any order of perturbation theory lead to the usual asymptotic states on 
the mass shell and no new states appear. From this point of view the theory remains 
unitary in physical space. The scattering amplitudes also behave safely decreasing with 
momenta, thus do not violating the unitarity requirement. 

The only source of conceptual problem is the presence of poles in the propagator in 
the complex momentum plane. They may appear in the Euclidean region for negative 
k"^ (as in D=5) or for complex fc^ (as in D=7). In the first case, they create a problem 
when integrating over Euclidean momenta and one might take the integrals in a sense 
of a principle value, while in the second case this does not happen, but in both the cases 
the presence of this type of poles might be interpreted as appearance of new (ghost) 
states. Note, however, that the form of a dressed propagator is not specific for the 
1/Nf expansion, but is typical of any perturbation theory. For instance, in the usual 
QED or any other renormalizable theory the dressed propagator has the same form 
and may obtain poles in the momentum plane (the Landau pole is an example). This is 
incompatible with the Kallen-Lehmann representation leading to a non-positive spectral 
function [17J or to a negative contribution from pole terms, which signals negative norm 
states exist. 

A usual way to avoid this problem is to consider the effective theory at low momenta 
below the Landau pole assuming that the situation is improved in a proper higher 
dimensional theory. Indeed, in the four dimensional QED the Landau pole exists at 
the energies above the Planck scale and one can safely use QED at lower momenta. In 
our case, however, due to power like behaviour it happens much earlier and severely 
constrains the region of validity of this theory. 

There are several attempts to build renormalizable effective quantum gravity using a 
kind of 1/iV expansion ^18j, where the role of an expansion parameter 1/A^ is played by 
the number space-time dimensions. The large D limit in this case is very similar to the 
large Nc planar diagram limit in the Yang- Mills theory considered by 't Hooft [12]. The 
technique similar to the 1/A^/ expansion is used also in [19], though no large parameter 
appears. The author, as in [20] , sums up the soft graviton corrections to the propagators 
of the scalar field to get an improved propagator which decreases faster than any power 
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of |/c^|. Though this partial resummation is similar to our expansion, the absence 
of an expansion parameter does not justify, to our mind, the selected set of diagrams. 
From this point of view the 1/N expansion is more consistent and contains the guiding 
line for such a selection. 
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